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Abstract - In this paper we study strange matter by investigating the stability window within 
the QMDD model at zero temperature and check that it can explain the very massive pulsar 
recently detected. We compare our results with the ones obtained within the MIT bag model and 
see that the QMDD model can explain larger masses, due to the stiffening of the equation of state. 
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Introduction. — Neutron star is a very dense object 
believed to be the remnant of a supernova explosion. Its 
estimated mass lies around (1 — 2M®), its radius is of the 
order 10 km and its temperature of the order of lO 11 ^ at 
birth, cooling rapidly to about 10 10 i^ by emitting neutri- 
nos. In most physical models for neutron stars, the star 
is composed only of hadrons and leptons, whose stellar 
matter is a mixture of degenerate neutrons, protons and 
electrons. In all stellar models, the structure matter de- 
pends on the assumed equation of state (EOS), which is 
derived from effective models. One of the uncertainties 
related to neutron star properties is the ground state of 
nuclear matter. In most models, hadrons are assumed to 
be the true ground state of the strong interaction. Since 
the proposal by Witten [Tj that strange matter (SM) may 
be the actual ground state of baryon matter at high densi- 
ties, many investigations have been conducted to verify the 
veracity of this hypothesis and the implications in several 
areas of physics, astrophysics and cosmology with possible 
consequences in the QCD phase diagram. 

SM was first considered in calculations obtained within 
the MIT bag model framework [2|. More sophisticated 
treatments for SM, based on the Nanbu-Jona-Lasinio [3] 
and the color flavor locked phase models [3l|4] also exist in 
the literature. 

In [5], a confinement mechanism was introduced by as- 
suming that the quark masses are density dependent. This 
model, named QMDD, was then applied to describe SM 
and the related quark star properties in [6]. In the very 
same paper, the authors pointed out that the results ob- 
tained with the QMDD model were quite different from 
the ones obtained with the MIT model. Subsequently, in 



[7J, the authors claimed that this difference was due to an 
incorrect thermodynamical treatment of the problem and 
recalculated the equation of state showing that an extra 
term is present in the energy density and pressure of the 
system. This extra term results from the dependence of 
the quark masses on the baryonic density. 

An important ingredient in the SM hypothesis is the sta- 
bility window, identified with the constant values of the 
model that are consistent with the fact that two-flavor 
quark matter must be unstable (i.e., its energy per baryon 
has to be larger than 930 MeV, which is the iron bind- 
ing energy) and SM (three-flavor quark matter) must be 
stable, i.e., its energy per baryon must be lower than 930 
MeV. It was also shown in [7J that the zero pressure den- 
sity does not correspond to the minimum of the energy 
per baryon, as is normally the case, because of the quark 
mass density dependence. In [8], the author has shown 
that the pressure at the density corresponding to the min- 
imum of the energy density can be zero if it is calculated in 
a self-consistent way along the energy density. That calcu- 
lation requires another extra term in the grand canonical 
thermodynamical potential, which mimics the MIT bag 
constant, but instead of being constant, it increases with 
the increase of the density. 

Still another calculation of the equation of state based 
on thermodynamics and also on the general ensemble the- 
ory was obtained in [9] . The authors claim that the extra 
term reported in [7J is correct in the expression of the 
pressure, but should not be present in the energy density 
equation. Similar kinds of density dependence are found 
in hadronic models [TOlITT] and used to describe protoneu- 
tron star properties in p~2 l fT3] . Specifically, in [TT][T3] , the 
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hadronic masses are density dependent while in [TU1TT2] the 
coupling constants of the model are density dependent. In 
both approaches, the pressure also carries a rearrangement 
term that appears due to the density dependence, while 
such a term cancels out in the energy density expressions. 
In [pj, due to quark confinement and asymptotic freedom, 
another prescription for the quark masses is used and the 
values for the quark current masses are somewhat differ- 
ent from the standart ones. Normally mo u = fnod = 5 
MeV and mos is of the order of 150 MeV. In [9], mo u = 5 
MeV, m od = 10 MeV and m 0s = 80 or 90 MeV. Depend- 
ing on the parametrization used, SM is completely stable 
or becomes metastable. 

More recently, the QMDD model was again revisited 
[14] . As already explained in this Introduction, the grand 
canonical thermodynamical potential is not only a func- 
tion of T, V and the chemical potential p, but also depends 
explicitly on the quark (or related baryon) density through 
its mass and this fact caused all the different treatments 
described below in the previous papers. As pointed out in 
[14] , the extra terms for the different treatments contradict 
each other and tackling thermodynamics self-consistently 
in the QMDD model seems to be a real problem. Another 
recipe is then proposed, considering an ideal quasi-particle 
system, where the mass can be both density and tempera- 
ture dependent, but is taken as constant at fixed values of 
To and pq. In this case, the standard ideal gas expressions 
are recovered and the extra terms do not appear. More- 
over, the pressure of the system is always positive and this 
poses a problem related to the stability of SM. In order to 
circunvent this difficulty, the authors have to consider the 
physical vaccum. 

In 2010, a measurement of the Shapiro delay in the radio 
pulsar PSR J1614-2230 yielded a mass of 1.97 ± 0.04M Q , 
the highest mass of a compact object ever measured. It 
is well known that the stellar mass that can be supported 
against gravitational collapse depends on the equation of 
state used to describe it. According to [16], very massive 
stars, with masses 2M & < M < 2.73M Q can be inter- 
preted as compact stars composed entirely of deconfined 
quarks and would be quark stars. 

The aim of the present work is to check whether this 
high mass pulsar can be described by the QMDD model. 
A necessary investigation that preceeds the use of the 
equation of state as input to the TOV equations [17] is 
the detailed consideration of the stability window. In the 
present paper we consider quarks with equal chemical po- 
tentials, closer to the situation in 56 Fe where matter is 
almost symmetric, than to the situation of charge neu- 
traliy usually analysed, where electrons, not present in 
nuclei, are included. Whenever possible, our results are 
compared with the MIT model. In face of all the dis- 
cussions on the thermodynamical consistency of different 
versions of the QMDD model, we restrict ourselves to the 
version obtained in [7] . 



Quark matter. — In this section we summarize the 
models we use to describe the properties of quark matter. 
The QMDD model is based on a phenomenological ap- 
proach [7] where the dynamical masses of the three light- 
est quarks scale inversely with the baryon number density 
[6], what mimics an interaction among the quarks: 
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where C is an energy density constant in the zero quark 
density limit, 
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is the baryon number density and m 0s s is the current mass 
the of the s quark. The energy density, the pressure and 
the quark density are respectively given by 
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where q = u,d,s, and m* is the effective quark 
mass. The distribution functions for quarks and 
anti-quarks are the Fermi-Dirac distributions f q ^ = 
1/ (1 + cxp [(£?* (p) =p fj, q ) /T]), with the chemical poten- 
tial for quarks (upper sign) and anti-quarks (lower sign) 

and E* (p) = ^Jp 2 + to* 2 . For T = 0, there are no an- 

tiparticles, the chemical potential is equal to the Fermi 
energy, and the distribution function for the particles is 

the usual step function: f qT = Q q (p 2 .- — p 2 ^j . 

In the description of compact stars, both charge neu- 
trality and /3 equilibrium conditions have to be imposed 
0: 



and 



2/O u = pd + p s + 3 (p e + Pu) 

fl s = p d = Pu + Pe, Pe = /V 



(6) 
(7) 



For the electron and muon pressure, energy density and 
densities we just replace q — > I in equations ([3]),(U]) and 
([5]), where I = e,p and 7/ = 2. 

We next summarize the main formulae at T = for 
both models used in this paper, the QMDD model [7] and 
MIT bag model 0. 
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The energy density, the pressure, the quark density and 
the baryonic density given in eqs. ©J® and ([5]) can be 
rewritten respectively as: 



p = H lj Br F ^ q )-B{c, m * q ) 
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where the dynamic term of confinement is given by: 
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Note that the term B (C, to*) arises from the derivatives 
of the grand canonical thermodynamical potential of a 
Fermi gas with respect to to* [16] . This density dependent 
term leads to a thermodynamical inconsistency, because 
the fundamental relation of thermodynamics £1 = —pV is 
violated. Although the model has this problem, we note 
that this density dependent term leads to an empirical 
model incorporation of the quark confinement at the low 
density regime and to asymptotic freedom at very high 
densities, as already emphasized in the Introduction. 

For the usual MIT bag model [2] the quark masses are 
fixed as m u = = 5 MeV, m s = 150 MeV and its 
expressions are dU),®, (fTO)) where to* is simply replaced 
by to and B (C, to*) by a bag constant B. 



Results and conclusions. — We now want to estab- 
lish the conditions under which SM is the true ground 
state [HE]- l n principle, the theory of the strong interac- 
tion should contain the answer to the question of whether 
strange matter is stable. Unfortunately, as we all know, 
QCD is still far from being completely solved. We need 
to rely on effective models to test the idea of SM. Hence, 
we next study strange matter via two different effective 
models and find a sizeable region of parameter space for 
which strange matter is stable. For the QMDD case, these 
parameters are C and the current strange quark mass TOo s 
(see equation ([1])). For the MIT bag model, these param- 
eters are the well-known B 1 / 4 and and current strange 
quark mass mo s . Strange matter is stable at zero pres- 
sure if its energy density per baryon density is lower than 

the energy per baryon of 56 Fe, =4] < 930MeV 

while the two- flavor quark matter must be unstable (i.e., 
its energy per baryon has to be larger than 930 MeV). 
Notice that, if finite size effects are taken into account, 
this value is around 4 MeV lower, i.e. 934 MeV [2], but 
the final stability window remains practically unchanged. 
The numbers arc discussed below. The energy per nucleon 
has been calculated numerically for 2QM and SM respec- 
tively. Therefore, SM is stable in the shaded region shown 
in Fig.Hl for the QMDD model and in Fig. [2] for the MIT 
model. The lower limit, vertical straight line, is due to the 
requirement that two flavor quark matter is not absolutely 
stable. In both figures one can see that we have allowed the 
strange quark mass to vary considerably. We have checked 
that the results arc slightly different if we consider matter 
with identical quark chemical potentials, corresponding to 
symmetric matter in the 2QM or charge neutral matter in 
/3-equilibrium, as expected in stars. For both models, the 
stable region is larger if the conditions given in eqs.© 
and (O are imposed. This is the situation which is usu- 
ally analysed in the literature, but one has to bear in mind 
that stable nuclear matter (as in iron) does not obey these 
conditions and no electrons are present. Actually, its pro- 
ton fraction is Y p = 0.46, very close to symmetric matter 
and this is the reason why we have chosen to analyse mat- 
ter with equal quark chemical potentials. If we consider 
finite size effects and take m,Q S = 150 MeV, the border of 
the stability line moves from C = 80.1 to 79.2 MeV/fm 3 
in the QMDD model and from B 1 / 4 = 155.1 to 155.8 in 
the MIT model. 

In Fig. [3^, we plot the equations of state (EoS) obtained 
for two constants in the MIT model and one situation in 
the QMDD model, chosen from the values which satisfy 
the SM stability condition. For these choices of param- 
eters, the QMDD EoS is reasonably stiffer. The MIT 
bag model presents a linear dependence of the pressure 
against the energy density as shown in Fig. |3j The pa- 
rameter B is just a linear parameter in the EOS for the 
MIT case. This is not true for the QMDD model, because 
the pressure as a function of the energy density is more 
complicated, as seen in Fig. |3J This non-trivial behav- 
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Fig. 2: SM stability window obtained with the MIT model at T = 0. The points in the figures show the maximum values for 
B 1/4 when m 0s = 150MeVis fixed. 
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ior is due to B (C, m*) in equation (fT2"]) . If we consider 
the thermodynamic treatment used in [5], where the zero 
pressure point corresponds to the minimum of the energy 
per baryon, the stability window modes to larger parame- 
ters, as seen in Fig. [3J3 and the equations of state become 
softer. 

As stated in the Introduction, in order to study hy- 
drostatic equilibrium of a family of quark stars we have 
to solve numerically the Tholman-Oppenheimer-Volkoff 
equations (TOV) |17j . The mass radius relation is ob- 
tained from the solution of these coupled differential equa- 
tions. 

The determination of small radii as the ones expected 
in pulsars from observations based on luminosity and tem- 
perature is full of uncertainties [18]. Hence, although X- 
ray astronomy has provided some simultaneous determi- 
nation of masses and radii from X-ray bursters [TS], they 
have to be taken with care. The integration of the TOV 
equations gives quark stars whose mass-radius relation is 
shown in Fig. 21 In this figure the horizontal lines are the 
masses of some well-known pulsars extracted from [151120] . 
The lower and upper limits of the masses and radii of 
EXO 0748-676 and 4U 1608-52 are also displayed. The 
shaded clouds refer to the lcr and 2er confidence ellipse of 
the results obatined in [19] for the EXO 1745-248. We 
notice that the QMDD model can certainly describe com- 
pact objects which are more massive than the MIT model, 
as already expected from its harder EoS. Nevertheless, it 
fails to describe pulsars with low radii, which can be de- 
scribed by the MIT model. One, however, has to bear in 
mind the uncertainties discussed in [18]. In Fig. [5] we 
show the maximum masses obtained in different points of 
the stability window for both models. They increase with 
the decrease of the strange quark mass and change vary 
little for the paramenters on the right border line. For 
the present EOS, the QMDD model yields sistematically 
larger maximum masses than the MIT model. 

All the equations of state obtained for the QMDD 
model discussed in the Introduction rely on derivatives 
of the grandcanonical thermodynamical potential, which 
are truncated at different arbitrary points. To obtain the 
correct expressions, an alternative calculation can be done 
via the energy- momentum tensor as [10|il2(ll3J . This work 
is in progress. 

In a forthcoming paper, we will analyse carefully the 
stability windown for protoquark stars described by quark 
matter at finite temperature (or fixed entropy). The choice 
of appropriate parameters in the search for stable SM at 
finite temperature systems requires, in addition to the in- 
vestigation presented in this paper, a carefull study of the 
free energy per baryon (J- = e — TS, where S is the en- 
tropy density of the system) , in contrast to previous works 
[2~iTf2"4] . where the energy density per baryon was used. Of 
course, the lower limits, given by the straight lines on the 
left of figures [T] and [5] remain the same because the 56 Fe 
ground state is a zero temperature system, where the bind- 
ing energy and the free energy per baryon coincide. 
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Fig. 3: In a) we show strange matter equations of state for 3 parameter choices and in b) we compare the stability window for 
two different versions of the QMDD model. 
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Fig. 4: Mass-radius relations for different strange stars. 




Fig. 5: Inside the flags we show the maximum masses (M/M Q ) for the corresponding parameter sets inside the stability window 
for a) the QMDD and b) the MIT models respectively. 
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